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Abstract 

We introduce an ad-hoc electrodynamics with advanced and retarded Lienard-Wiechert interac- 
tions plus the dissipative Lorentz-Dirac self-interaction force. We study the covariant dynamical 
system of the electromagnetic two-body problem, i.e., the hydrogen atom. We perform the linear 
stability analysis of circular orbits for oscillations perpendicular to the orbital plane. In particular 
we study the normal modes of the linearized dynamics that have an arbitrarily large imaginary 
eigenvalue. These large eigenvalues are fast frequencies that introduce a fast (stiff) timescale into 
the dynamics. As an application, we study the phenomenon of resonant dissipation, i.e., a motion 
where both particles recoil together in a drifting circular orbit (a bound state), while the atom 
dissipates center-of-mass energy only. This balancing of the stiff dynamics is established by the ex- 
istence of a quartic resonant constant that locks the dynamics to the neighborhood of the recoiling 
circular orbit. The resonance condition quantizes the angular momenta in reasonable agreement 
with the Bohr atom. The principal result is that the emission lines of quantum electrodynamics 
(QED) agree with the prediction of our resonance condition within one percent average deviation. 
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I. INTRODUCTION 



In this paper we experiment with the stability analysis of the circular orbits of the electro- 
magnetic two-body problem. The main motivation is to understand the complex dynamics 
described by the electromagnetic equations of motion, that involve delay and third deriva- 
tives. We give a method to derive the linearized equations of motion in the neighborhood 
of the circular orbits of this implicitly-defined dynamical system with delay. We introduce 
an ad-hoc electromagnetic-like setting that uses advanced and retarded Lienard-Wiechert 
interactions plus the dissipative Lorentz-Dirac self-interaction force [1], henceforth called 
the dissipative Fokker setting (DFS). We study in detail a specific feature of the tangent 
dynamics of the circular orbits of the two-body problem; The stiff normal modes of the 
linearized dynamics, that have an arbitrarily large imaginary eigenvalue. Last, we discuss 
an application for the hydrogen atom and the surprising predictions of stability analysis 
within the DFS; We predict several features of the Bohr atom [2] with high precision and 
qualitative detail. A subset of the emission lines predicted by the DFS agrees with the lines 
of quantum electrodynamics (QED) within one percent average deviation. There is also a 
surprising body of qualitative agreement with QED; (i) the emitted frequency is different 
from the orbital frequency, (ii) the stable orbits of the DFS have angular momenta that are 
multiples of a basic angular momentum. This basic angular momentum of the DFS agrees 
well with Planck's constant and depends only logarithmically on the mass of the heavier 
particle. 

Dirac's 1938 fundamental work [1] on the electrodynamics of point charges gave complex 
and stiff delay equations that were seldom studied. Among the few models studied within 
Dirac's theory, the system of Eliezer's theorem [3-6] revealed a surprising dynamics; An 
electron moving in a Coulomb field with inclusion of self-interaction can never fall into the 
center of force by radiating energy. The result was generalized to motions in arbitrary attrac- 
tive potentials [4], as well as to tridimensional motions with self-interaction in a Coulomb 
field [5, 6], finding that only scattering states are possible. Since our model has Eliezer's 
problem as the infinite-mass limit, a finite mass for the proton is essential for a physically 
meaningful dynamics; If the proton has a finite mass, there is no inertial frame where it 
rests at all times, and this in turn causes delay because of the finite speed of light. It is 
widely known that QED gives a satisfactory and precise description of atomic physics, but 
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the same is not popularly thought about atomic models based on classical electrodynamics. 
Since dynamical studies are still missing, clearly this complex dynamics needs to be investi- 
gated beyond our preliminary findings. Even though we are not trying to replace QED, our 
understanding of this two-body dynamics might prove useful for atomic physics and perhaps 
we can understand QED as the effective theory of this complex stiff dynamics with delay. 
We shall describe the two-body motion in terms of the familiar center- of-mass coordinates 
and coordinates of relative separation, defined as the familiar coordinate-transformation that 
maps the two-body Kepler problem into the one-body problem with a reduced mass. We 
stress that in the present relativistic motion the Cartesian center-of-mass vector is not ig- 
norable, and it represents three extra coupled degrees-of-freedom. We introduce the concept 
of resonant dissipation to exploit this coupling and the many solutions that a delay equation 
can have. Resonant dissipation is the condition that both particles decelerate together, i.e., 
the center-of-mass vector decelerates, while the coordinates of relative separation perform an 
almost-circular orbit, despite of the energy losses of the metastable center-of-mass dynamics. 

Historically Nordstrom [7, 8] suggested the use of advanced and retarded potentials in 
atomic physics already in 1920, but the self-interaction theory was problematic in 1920 and 
the idea disappeared. The theory of nonlinear dynamics was not out yet in 1938 when 
Dirac's theory for the electrodynamics of point charges appeared [1], neither in the Glorious 
days of the twentieth century physics [9], such that the our present experiment is a new 
application of modern nonlinear dynamics. Advanced interactions appeared again in 1945, 
when Wheeler and Feynman [10, 11] gave an electrodynamics based on the postulate that 
every field is produced by charges located somewhere [12]. The theory was called action- 
at-a-distance electrodynamics [10, 11], a theory where the isolated two-body problem is 
defined by Fokker's action 

S F = - J m 1 ds 1 - J m 2 ds 2 -e 1 e 2 J J S(\ |xi - x 2 | | 2 )xi • ± 2 ds 1 ds 2 , (1) 

with Xj , Si, nii and e, representing the four-position, the proper time, the mass and the 
charge of particles % — 1,2 respectively. In Eq. (1) the dot indicates the Minkowski scalar 
product of four-vectors and double bars stand for the four- vector modulus [10, 11]. Due 
to the similarities with the equations of motion of the DFS, the dynamical studies of the 
action-at-a-distance theory are relevant for the present work. For example, in the collision 
of two electrons with equations of motion determined by Eq. (1), the solution is determined 
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by initial position and velocity only, as proved in Ref [13] (a Banach-to-Banach contraction 
mapping proof for nonrunaway orbits). This suggests that we are dealing with a perfectly 
causal and well-posed dynamical system dressed in unusual form [14]. Driver's result [13] 
suggests that a dynamics with advance and delay is well-posed in the same way. The DFS 
presents exactly the same neutral-delay mathematical problem of any electromagnetic-like 
model, as for example the problem with retarded-only fields of Refs.[15, 16]. Fokker's action 
of Eq. (1) is used here to derive the sector of the DFS equations of motion determined 
by the semi-sum of Lienard-Wiechert fields. Last, advanced interactions appeared again in 
another work of Eliezer; a generalization of Dirac's covariant subtraction of electromagnetic 
infinities [17]. The resulting generalized electromagnetic settings include the advanced inter- 
actions naturally, and provide a testbed for future studies in electrodynamics [17]. Here we 
shall keep to the DFS as a generic electromagnetic-like example. 

The road map for this paper is as follows; In Section IV we give the main technical 
part of the paper; We outline an economical method to derive the tangent dynamics of the 
circular orbit based on a quadratic expansion of the implicit light-cone condition. In this 
Section we also take the stiff limit of the linear modes of the tangent dynamics. In Section 
V we give an application to atomic physics, by discussing a necessary condition for the 
state of resonant dissipation; This condition is heuristically expressed by a simple resonance 
condition that predicts the correct atomic scales. The earlier sections are a prelude to 
Section IV. In Appendix A we discuss how the DFS can be fit into Dirac's electrodynamics 
of point charges. Section II is a review of the circular orbit solution, to be used in Section 
IV and in Section V. In Section III we build familiarity with Fokker's action of the action- 
at-a-distance electrodynamics as a prelude to the quadratic expansions needed for the linear 
stability analysis of Section IV. Last, in Appendix B we discuss the soft normal modes of 
the tangent dynamics and in Section VI we put the conclusions and discussion. 



II. THE CIRCULAR ORBIT SOLUTION 



In this Section we review the circular-orbit solution of the isolated electromagnetic two- 
body problem of the action-at-a-distance electrodynamics [18, 19], to be used as the un- 
perturbed orbit. For the isolated electromagnetic two-body problem, the tangent dynamics 
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studied in the next section is straight Lyapunov stability analysis. In the DFS there is also a 
very small force along the orbital plane of the circular orbit, such that a non-drifting circular 
orbit is not a solution of the equations of motion. In the DFS the tangent dynamics is the 
starting point of a perturbation scheme to impose that the resulting dynamics is a drifting 
circular orbit (the state of resonant dissipation). 

We use the index % = 1 for the electron and i = 2 for the proton, with masses mi and 
m 2 respectively, as of Eq. (1). We henceforth use units where the speed of light is c = 1 
and ei = — e 2 = — 1 (the electronic charge). The circular orbit is illustrated in Fig. 1; A 
motion of the two particles in concentric circles with the same constant angular speed and 
along a diameter. This dynamics satisfies the time-symmetric problem of Fokker's action 
(1) because the symmetric contributions from future and past generate a resulting force 
normal to the velocity of each particle [18, 19]. The details of this relativistic orbit will 
be given now; The constant angular velocity is indicated by Q, the distance between the 
particles in light-cone is r& and 9 = Qr& is the angle that one particle turns while the light 
emanating from the other particle reaches it (the light-cone time lag). The angle 9 is the 
natural independent parameter of this relativistic problem. Each particle travels a circular 
orbit with radius and scalar velocity defined by 

n = b in , (2) 

r 2 = b 2 r b , 

and 

Vl = Q n = Oh, (3) 
v 2 = fir 2 = 9b 2: 

for the electron and for the proton, respectively. The condition that the other particle turns 
an angle 9 during the light-cone time lag is [18] 

&? + &i + 2&i&2Cos(0) = l, (4) 

and is henceforth called the unperturbed light-cone condition. In Appendix B we calculate 
b\ and b 2 in a power series of 9 up to the fourth order. Last, because of the rotational 
invariance of Fokker's action, there is a conserved angular momentum perpendicular to the 
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plane of the orbit, that is evaluated in Ref. [18] to be 

1 + v x v 2 cos(6>) 

lz ~ + Vl v 2 sm(ey [ ) 

where the units of l z are e 2 /c, just that we are using a unit system where e 2 = c = 1 
[18]. Equation (4) restricts b\ and b 2 to be less than one such that for small values of 6 
the angular momentum of Eq. (5) is of the order of l z ~ 9" 1 . For orbits in the atomic 
magnitude, l z ~ O^ 1 is about one over the fine-structure constant, a^ 1 = 137.036. It is 
curious to notice that each (advance/retarded) interaction term of Fokker's action, Eq. (9 ), 
evaluates exactly to \Vtl z along a circular orbit, with l z given by Eq. (5). This combination 
of angular momentum times the orbital frequency is reminiscent of the formal maneuvers of 
quantum mechanics. 



III. FOKKER'S ACTION 



We use Fokker's action in this work as a means to derive the sector of the DFS equations of 
motion determined by the semi-sum of the Lienard-Wiechert potentials. In the following we 
discuss the Lagrangian formalism of Fokker's action (1) as an introduction to our economical 
method to obtain the tangent dynamics by expanding this action to quadratic order. The 
delta-function of Fokker's action (1) contains the retarded and the advanced light-cone 
contributions, and it is convenient to separate those two parts by factoring the argument of 
the delta function as 

(ti - h) 2 - r\ 2 = [h -t 2 - r 12 ] [h -t 2 + r 12 ], (6) 

where r 12 stands for the Cartesian distance between particle 1 at time t\ and particle 2 at 
time t 2 and each factor on the right-hand side of Eq. (6) is related to the advanced and the 
retarded light-cones of particle 1 respectively. The delta-function of a product argument 
is a sum of two delta-functions each multiplied by the respective Jacobian, such that the 
interaction term of Fokker's action (1) can be written as 

VA = [ -^-Sih - t 2 - r 12 )(l - v ± ■ v 2 )dt 1 dt 2 (7) 
J 2r i2 

+ / ~ — 8(t 1 -t2 + r 1 2)(l-vi-v 2 )dtidt2, 
J 12 
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where vi henceforth stands for the Cartesian velocity of particle 1 at time t\ and v 2 hence- 
forth stands for the Cartesian velocity of particle 2 at time t 2 . We henceforth use the dot to 
indicate the scalar product of two Cartesian vectors, as already used in Eq. (7). Integration 
of each term of Eq. (7) over t 2 brings out another Jacobian factor and yields 

VA= [ 1 (l-^-^) dti + I \ dtu ( 8 ) 

J 2r 12 (l + n 12a • v 2a ) J 2r 12 (l - n 126 • v 26 ) 

where n 12a is a unit vector connecting the advanced position of particle 2 at time t 2 to the 
position of particle 1 at time t±, vector n^f, is a unit vector connecting the retarded position 
of particle 2 at time t 2 to the position of particle 1 at time tiand v 2a and v 2 & stand for the 
velocity of particle 2 at the advanced and retarded time i 2 respectively. Equation (8) is the 
most useful form of Fokker's interaction for our purposes. Notice that each term of Eq. (8) 
can be cast in the form 

1 f d-v.-v,) -1 I _ vi _ 



2 J r 12 (l + ilia) 1 2 

where V and A are the Lienard-Wiechert scalar potential and the Lienard-Wiechert vector 
potential respectively. We have introduced the quantity c = ±1 in the denominator of Eq. 
(9) such that c = 1 represents the advanced interaction while c = — 1 represents the retarded 
interaction. The quantities of particle 2 in Eq. (9) are to be evaluated at a time t 2 defined 
implicitly by 

t 2 = t 1 + — , (10) 

c 

where c = ±1 describes the advanced and retarded light cones, respectively. Because of this 
decomposition of Fokker's interaction into V and A parts, we henceforth call Eq. (9) the 
VA interaction. A derivation of the Lienard-Wiechert potentials from Fokker's action and 
details such as the Darwin approximation are found in Ref.[20]. 

The stiff limit is determined by the largest-order derivative appearing in the linearized 
equations of motion of Appendix A. In this approximation, the contribution of the self- 
interaction force to the linearized dynamics about a circular orbit is simply given by the 
Abraham-Lorentz -Dirac force 

Frad = jja. (11) 

The contribution of the other smaller terms will be given elsewhere. 
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IV. LINEAR STABILITY ANALYSIS 



In this Section we study the linear stability analysis of the circular orbits for displacements 
perpendicular to the orbital plane, henceforth called the z-direction. We give an economical 
method to obtain these equations of tangent dynamics by expanding the implicit light-cone 
condition up to quadratic order. We start from the equations of motion of the isolated 
system, which are derived from Fokker's action (1) and yield the Lienard-Wiechert fields in 
the half-retarded plus half-advanced combination. This linearized z dynamics is uncoupled 
from the planar dynamics, and the linearized equations can de derived without the use of a 
symbolic software, as we explain in the following. The Cartesian coordinates of a transversely 
perturbed circular orbit are defined by 

Xk + Wk = r b d k exp(iQt), (12) 
Xk-iVk = r b d* k ex.p(-ittt), 
Zk = r b CSZ k , 

where k = 1 for the electron and k = 2 for the proton, Z k is the small transverse perturba- 
tion, d\ = bi and d 2 = —b 2 are defined from the two real parameters of Eq. (2) and fl is the 
orbital frequency defined above Eq. (2). Last, in Eq. (12) C and S are defined by 

C = l + &!& 2 2 cos(0), (13) 

and 

S= l + 6iMsin(0). (14) 

We henceforth introduce a scaled time r = fit. The linear stability analysis involves 
expanding the equations of motion to linear order in Z k , which in turn is determined 
by the quadratic expansion of Fokker's action in Z k . The main tool for expanding this 
quadratic action is the perturbed light-cone condition, Eq. (10), about the circular orbit 
(where ryi = r b is the constant circular lag). We introduce a function (p of the Z\ and Z 2 
perturbations by expanding the light-cone time t 2 as 

+ ? + & (15) 

In the following we calculate this homogeneous functional (p of Z\ and Z 2 up to quadratic 
order. The distance r 12 entering Eq. (10) is to be evaluated from the position of particle 1 
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at time t±, to the position of particle 2 at the time t 2 defined implicitly by Eq. (15). The 
coordinates of particle 2 at the time t 2 are defined implicitly by 

%2 + iy2 = r b d 2 exp(iri + ic9 + iip), (16) 
%2 — W2 = rbdl exp(— vr\ — ic9 — iip), 

z 2 = r b CSZ 2 {r 1 + c9 + ip) = r b CSZ 2c , 

where c = 1 for the advanced time t 2 and c = — 1 for the retarded time t 2 . Notice that Eq. 
(16) defines the coordinates implicitly, because tp is a function of the deviations Z\ and Z 2 . 
Even though 9 is small in applications of atomic physics, we stress that one should never 
expand in powers of 9; the correct infinitesimal quantity of the linear stability analysis is the 
size of the deviations from circularity and their homogeneous functions such as p (expanding 
in 9 produces the Darwin approximation [20]). This non-analyticity will become clear after 
we show that the logarithm of 9 appears. We therefore expand the advanced/retarded 
position Z 2c of particle 2 at the scaled time T\ + c9 + <p in a Taylor series in <p about the 
advanced/retarded position t 1 + c9. It turns out that only the zeroth-order term appears 
in the action up to quadratic order. Because of this, the linearized equations involve only a 
constant shift, a considerable simplification. Substituting t 2 of Eq. (15) together with the 
position (16) of particle 2 into Eq. (10) and using Pythagoras theorem for the distance r 12 
from particle 1 at time ti to particle 2 at time t 2 yields 

r\ 2 = (n + n^-y = rl + rlC 2 S\Z 1 - Z 2c f. (17) 

Notice that the Z variations decouple from the planar variations because there is no mixed 
linear term of Z times a linear perturbation of the planar coordinate in Eq. (17); These are 
naturally separated by Pythagoras theorem. The planar perturbations enter in Eq. (17) as 
an added quadratic form, as given in the next section. It is convenient to define another 
function $ by ip = 9cCS§, such that Eq. (17) is a quadratic equation of $ and the regular 
solution up to second order in Z\ and Z 2c is 

C 9 

$> = —{Z 1 -Z 2c f. (18) 

The coordinate Z 2 appears evaluated at the advanced/retarded time in Eq. (18), and to 
obtain the action up to quadratic terms it is sufficient to keep the first term Z 2c = Z 2 {j\ + 
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c9 + ip) ~ Z 2 (ti + c9). Using the 2— perturbed orbit defined by Eq. (16) to calculate the 
numerator of the VA interaction of Eq. (9) yields 

(1 - Vl • v 2c ) = 1 + 9 2 cos(0)&i&2 - 9 2 C 2 S 2 Z 1 Z 2c , (19) 

and the denominator of the VA interaction of Eq. (9) is 

7-12(1 + n 12c • v 2c /c) = r b S[l + C$ + 9cC 2 S(Z 1 - Z 2c )Z 2c \. (20) 

Notice that the quadratic term Z 2c Z 2c on the right-hand side of Eq. (20) can be dropped 
because it represents an exact Gauge that does not affect the Euler-Lagrange equations of 
motion, such that 

r 12 (l + n 12c • v 2c /c) » r h S[l + C$ + flcC 2 ^^], (21) 

where the equivalence sign pa henceforth means equivalent up to a Gauge term of second 
order. Even if a quadratic Gauge term appears in the denominator, in an expansion up to 
quadratic order it would still produce a Gauge and therefore it can be dropped directly from 
the denominator. One should be careful not to do this with linear Gauges, which appear 
only in the planar stability analysis to be considered elsewhere. In this way, the expansion 
up to second order of the VA interaction of Eq. (9) is simply 

VA pa (^){1 - 9 2 CS 2 Z 1 Z 2c - < ^-{Z l - Z 2c ) 2 - 9cC 2 SZ 1 Z 2c }. (22) 

Last, we need the kinetic energy along the z-perturbed circular orbit, which we express in 
terms of Z\ of definition (12) as 



Tx = -m iy Jl - v 2 = ~yjl - TiCWZ 2 , (23) 

where the dot means derivative with respect to the scaled time r, 7^ = a/1 — v \ , and we 
have used Vtr b = 9. The expansion of Eq. (23) up to second order is 

r '=te ){ ^r + 2 z ' + ->- < 24 > 

where e\ = mir b 'j 1 9 2 CS 3 is calculated with Eq. (55) to be 

ei = ^{[C 2 + 9 2 S(S - l)](6i + 6 2 cos(^)) +S(9sm(9) - 9 2 cos(9))b 2 } . (25) 
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We are ready to derive the Euler-Lagrange equation of motion for particle 1 of the isolated 
two-body problem using the quadratic Lagrangian 



L 1 = T 1 + VA C=1 + VA c= -l (26) 



This equation of motion is 



/nr2 q f)C^ S* 9 2 C i S' 2 

e x Z x = -^(2^i - Z 2+ - Z 2 _) - ^(Z 2+ - Z 2 _) - ~^(Z 2+ + Z 2 _). (27) 

Notice that the term on the left-hand side of Eq. (27) can be written as 

e 1 Z 1 = r^m^CSZ, = r 2 b S 2 ^, (28) 

which is proportional to the force along the z-direction. According to the prescription of 
the DFS, we shall add the following self-interaction term to the right-hand side of Eq. (27) 

r 2 b S 2 F rad = ^CS 3 Zu (29) 

where the triple dot means three derivatives with respect to the scaled time and we have 
used Eq. (11). The full linearized equation of motion for Z 1 is 

e x Z x = -CS r z\ - ^(2^i - Z 2+ - Z 2 _) - ~^(Z 2+ - Z 2 _) - ~^{Z 2+ + Z 2 _). (30) 

The linearized equation for Z 2 is completely analogous and is obtained by interchanging 
Z\ by Z 2 and ei by e 2 in Eq. (30). The general solution of a linear delay equation can 
be obtained by Laplace transform [22] and is a linear combination of the following normal 
mode solutions. A normal mode solution is obtained by substituting Z\ = Aexp(pr) and 
Z 2 = B exp(pr) into the two linearized equations, and requires the vanishing of the following 
2x2 determinant 



detZ = 



(31) 



c 2 s + e lP 2 - lcs 3 e 3 P 3 g(o, p ) 

G(9,p) C 2 S + e 2 p 2 -ICS 3 9 3 p 3 

where G(9,p) = (C 2 S — C S 2 9 2 p 2 ) cosh(p9) + C 2 Sp9smh(p9). Two kinds of limits are in- 
teresting for the infinite-dimensional formal collection of normal modes of Eq. (31); (i) the 
four soft Coulombian modes obtained by expanding Eq. (31) in powers of 9 for small values 
of p, as discussed in Appendix B, and (ii) the stiff limit obtained when p9 is large, such that 
the hyperbolic functions of the G(9,p) acquire a large magnitude [23]. In the following we 

11 



use the zeroth-order term of the expansion for b\ and 62 given in Appendix B to evaluate 
the determinant (31) : 



ci = — + 0(9 2 ), (32) 
m 2 

«, = ^ + 0(«"). 
mi 

C = l + 0(# 2 ), 
S = l + O{0 2 ). 

For small 9, the second-order and higher even-order terms of Eq. (32) give only a small 
correction. Substituting Eq. (32) into Eq. (31) and defining p = X/9, we obtain 

-^(det Z) = l- 2 -B 2 \ + ^-£Va 2 - ^[(1 - A) c °sM A ) " T sinh(A)] 2 , (33) 
MA 3 9M M A A y n y 1 

where we have dropped small 0(9 2 ) terms. The stiff-mode condition defined by Eq. (31) 
(det Z = 0) is 



1 - 2 f 2 \ + ^ 4 A 2 - ^[(1 - ^ + ^) cosh 2 (A) + i(l - 1) sinh(2A)] = 0. (34) 



For future reference we give also the stiff limit for the ^-tangent dynamics without the 
self-interaction terms, which is obtained from Eq. (27) and the corresponding equation for 
particle 2 

1 - ^[(1 - -L + ^) cosh 2 (A) + 1(1 - 1) sinh(2A)] = 0. (35) 
V. THE STIFF STABILITY OF THE HYDROGEN ATOM 



We are interested in finding motions where the particles recoil together while staying in 
the neighborhood of a drifting circular orbit, i.e., the state of resonant dissipation. The 
need for a resonance becomes obvious in the following perturbative scheme; (i) We take 
the circular orbit as the unperturbed state, (ii) We substitute the circular orbit plus a 
perturbation into the equations of motion of the DFS and take the linearized equations of 
motion. The circular orbit is not an exact solution of the DFS equations of motion, because 
of the small forcing coming from the third derivatives. This perturbative scheme yields linear 
delay equations with a small forcing term along the orbital plane. It is then possible to show 
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by averaging [26] that a weakly-accelerated drifting circular orbit is never a solution to these 
linear equations. Therefore, a bifurcation of the circular orbit must happen and a nonlinear 
term must be important to balance the small dissipative forcing, if the state of resonant 
dissipation is to be attained. In the following we postulate that this resonance happens at 
a quartic order. By inspection, one finds that only resonance conditions involving the stiff 
modes can be satisfied in the atomic magnitude. In the following we study the consequences 
that along some special circular orbits such balancing mechanism is established by the 
existence of a quartic resonant constant of motion. To discuss this stability by resonance 
we need some results of the tangent dynamics along the orbital plane. This more elaborate 
tangent dynamics is derived in a way analogous to Section V and shall be given elsewhere, 
here we give only the main results. The stiff-limit for the equal-mass two-body problem 
with retarded and advanced fields is studied in Ref. [23], and in the following we give the 
generalization of these results for the case of arbitrary masses. 

Up to linear order, the tangent dynamics along the orbital plane is decoupled from the 
^-dynamics of Section V. To study this planar tangent dynamics, it is convenient to describe 
the orbit along the z = plane using gyroscopic coordinates 



where r] k and £ fc are complex numbers defining the perturbation of the circularity and the 
dk are defined below Eq. (12). Because x k and y& are real, we should have r\ k = £* k but a 
convenient way to minimize the quadratic functional of Fokker's action is to treat i] k and £ fc 
as independent functions. To fix ideas we start from the stability of the isolated two-body 
system, and again we define the normal- mode eigenvalue by XQ/9, i.e., every coordinate 
perturbation oscillates in time as exp(\Qt/9) ( A is an arbitrary complex number). The 
limiting form of the planar characteristic equation for the isolated different-mass case is 



where \i is the reduced mass and M = mi + m 2 (for the equal-mass case, our general 
Eq. (37) reduces to Eq. (15) of Ref. [23] ). Along circular orbits both the planar and 
the perpendicular linearized equations share the same limiting characteristic Eq. (37), as 
can be checked with Eq. (34). For hydrogen (fi/M) is a small factor of about (1/1824). 



%k + iyk = r b exp(iQt) [d k + rj k ], 
x k - iyk = r b exp(-iQt)[d k + f fc ], 



(36) 




(37) 
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It is important to understand the structure of the roots of Eq. (37) in the complex A 
plane, specially for 6 of the order of the fine structure constant. The very small parameter 
^ ~ 1CT 13 multiplying the squared hyperbolic cosine on the left-hand side of Eq. (37) 
determines that a = | Re(A)| ~ In(^Jj^). For the first 13 excited states of hydrogen this 
a is in the interval 14.2 < \a\ < 18.2. The imaginary part of A can be an arbitrarily large 
multiple of 7r, such that the general solution to Eq. (37) is 

A = ±(<7 + mq), (38) 

where q is an arbitrary integer. The plus or minus sign of Eq. (38) is related to the 
time-reversibility of the isolated two-body system, a symmetry that is broken by radiation. 
This same exact phenomenon happens for the z direction. Next we include the dissipation 
of the DFS, i.e., the Lorentz-Dirac self-interaction, a calculation performed by adding the 
self-interaction force to the equations of motion of the isolated system. Here we give only 
the characteristic planar equation up to 0(1/ A 4 ) 

(1 + ^ + ^)(^) cosh 2 (A) = 1 - Va + Va 2 + (I + sinh(2A) + ... (39) 

It is remarkable that Eqs. (39) and (34) differ only at the terms of 0(1/A) and at the terms 
of type 6* 4 A 2 , which describe small corrections for a in the atomic range. The linear term on 
the right-hand side of Eqs. (34) and (39) with the 2/3 coefficient is due to the self-interaction 
force. This dissipative term breaks the time-reversal symmetry of Eq. (37), and the roots 
of Eqs. (34) and (39) no longer come in plus or minus pairs. Let X xy be a root of Eq. (39) 
with positive real part and \ z be a root of Eq. (34) with a negative real part. In the stiff 
limit these are both near one of the limiting roots (38) and can be expressed as 

K y = (cr + -nqi + itx), (40) 
\ z = — (cr + nqi + ^2), 

where the small perturbations ei and e 2 are so far two arbitrary complex numbers. The 
second order balancing process studied here involves the interaction of a z mode with a 
planar mode, in the same way used in Refs. [15, 16]. This is because if the atom is to 
recoil like a rigid body, one expects the fast dynamics to encircle the circular orbit with fast 
spinning motions of balanced amplitude. 
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We henceforth assume heuristically that the state of resonant dissipation is formed in a 
bifurcation involving perturbations along two special linear modes of the tangent dynamics. 
We take a perpendicular normal mode of Eq. (34) and a planar normal mode of Eq. (39), 
with eigenvalues X z and X xy respectively. The coordinate of the planar normal mode is 
a linear combination of the four rj^ gyroscopic coordinates: u = a lk r\ k + &ifc£ fc , while the 
coordinate of the perpendicular z normal mode is Z = b\Z\ + b^Z2- Using the normal mode 
conditions 9ii = Q\ xy u and 9Z = Q\ Z Z one can show that the quadratic form uZ is a 
complex amplitude that oscillates harmonically with the beat frequency (X xy + \ z )Vt/9 = 
i{e\ — e2)Q/9. Our resonance condition is to choose these two eigenvalues such that 

Re(\ xy + X z ) = 0. (41) 

Condition (41) avoids that the modulus of the amplitude uZ has an exponential growth. We 
shall see that condition (41) is satisfied only for special discrete values of 9. Since condition 
(41) must be satisfied, we henceforth assume that t\ and €2 are real numbers, as any excess 
real part in Eq. (40) can be absorbed in the definition of a. Condition (41) is also the 
necessary condition to construct a resonant constant in the neighborhood of the circular 
orbit; Because Fokker's action is real, A* and \* xy are also eigenvalues to Eqs.(34) and (39) 
respectively, with complex conjugate normal mode coordinates. Condition (41) then implies 
the usual necessary condition for a resonant constant 

Ky + K + Ky + K = 0, (42) 

as discussed in Refs. [16, 24]. Using these complex conjugate normal-mode coordinates and 
Eq. (41), one can show that the following quartic form is a constant of the motion up to 
higher order terms [16, 24]: 

C =\u\ 2 \Z\ 2 + ... (43) 

The quartic function of Eq. (43) is constant because it is the squared modulus of the 
harmonic amplitude uZ = v / Cexp(i(ei — e 2 )Vlt/9). This necessary condition and the con- 
tinuation of the leading term (43) to an asymptotic series is discussed in Ref. [16]. 

The root- searching problem of Eq. (41) is well posed and for each integer q conditions 
(34) and (39) together with Eq. (40) determine a unique 9 as a function of q , i.e., 9 is 
quantized by the integer q that appears naturally in Eq. (40). An asymptotic solution to 
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condition (40) can be obtained by expanding Eqs. (34) and (39) up to quadratic order in 
ei and €2 while treating a as an approximate constant. This approximation determines the 
following discrete values for 9 



6(7rV - a 2 ) 
a{7i 2 q 2 + a 2 ) 2 ' 



(44) 



and 



(ei - e 2 ) 



4Trq(3a 2 - TX 2 q 2 ) 
a{a 2 + n 2 q 2 ) 2 



(45) 



According to QED, the circular Bohr orbits have maximal angular momenta for that quan- 
tum number and a radiative selection rule ( Al = ±h) restricts the decay from level k + 1 
to level k only, i.e. circular orbits emit the first line of each spectroscopic series (Lyman, 
Balmer, Ritz-Paschen, Brackett, etc.), henceforth called the QED circular line. We have 
solved Eqs. (34), (39) and (40) with a Newton method in the complex A plane. Every angu- 
lar momentum l z = Q~ x determined by Eq. (41) has a value in the correct atomic magnitude 
( d^ 1 > 137.0 ); The first resonance appears at q = 5 for Q~ x = 252 A and the minimum value 
d^ 1 = 48.52 is attained at q — 7, then # _1 increases monotonically with q. The subset of 
Table 1 has frequencies wdf surprisingly close to the QED lines. These lines are for q approx- 
imately equal to an integer multiple of the integer part of 2a. We conjecture here that among 
the resonances satisfying the necessary condition (41), only some have \u\ 2 depending on the 
translation-invariant quantities (£ x — £ 2 ) an d ~~ V2) to allow a recoiling translation [26]. 
In our description the emission mechanism is at a frequency equal to the orbital frequency Q 
corrected by the frequency of the complex amplitude uZ defined above Eq. (43), as we ex- 
plain below. The numerically calculated angular momenta l z = for this select subset are 
given in Table 1, along with the orbital frequency in atomic units (137 3 f2)//i = (1376*) 3 , the 
QED first frequency of the series in atomic units wqed = |(p — (FPTp)> anc ^ ^ ne f rec L uenc y 
predicted by the dissipative Fokker model wdf = (1370) 3 + 137 3 # 2 (e! — e 2 ). We list only the 
first 13 lines, which are the experimentally observable, but we tested the agreement of the 
numerical calculations of the Newton method with up to the 40 th circular line predicted by 
QED. Beyond that, the asymptotic formula (44) shows that the agreement is essentially for 
any integer k because substitution of q — [2a]k into Eq. (44) yields 




(46) 
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to be compared with the 137.036 of QED. The agreement for any integer k suggests that 
Eqs. (34) and (39) describe a linear problem that is equivalent to Schroedinger's equation 
(linear operators with the same spectrum are equivalent). 



l z = 


(137#J 


Wqed 


w D F 


Q 


161.94 


6.054x10 


o ^7 r A , , 1 /"I — 1 

3.750x10 


3.655x10 


32 


o o o r o 

283.52 


1 1 OO ^ , 1 O — I 

1.128x10 


6.944x10 


6.774x10 


r r 

55 


398.06 


A O TT ^ , 1 O — 2 

4.077x10 


2.430x10 


o a o o v * i o — 2 

2.462x10 


76 


r o ( \ oo 

520.29 


1.826x10 


i 1 o r ^ , i o — 2 

1.125x10 


1.110x10 


OO 

98 


638.53 


9.876x10 


6.111x10 


OO O ^ , 1 O — 3 

6.038x10 


1 1 O 

119 


T r o o 'i 

752.27 


6.039x10 


o n o r - , 1 o — 3 

3.685x10 


3.710x10 


1 OO 

139 


872.68 


3.868x10 


2.406x10 


2.387x10 


160 


988.16 


2.664xl0- 3 


1.640xl0~ 3 


1.650xl0- 3 


180 


1110.15 


1.879xl0~ 3 


1.173xl0~ 3 


1.168x10-3 


201 


1226.95 


1.392xl0- 3 


8.678X10" 4 


8.677X10" 4 


221 


1344.30 


1.058xl0- 3 


6.600xl0~ 4 


6.615xl0" 4 


241 


1462.14 


8.226xl0~ 4 


5.136X10" 4 


5.153xl0" 4 


261 


1580.44 


6.513xl0" 4 


4.076xHT 4 


4.090xHT 4 


281 



Table 1: Numerically calculated angular momenta l z = 9~ x in units of e 2 /c, the orbital 
frequencies in atomic units (1376*) 3 , the circular lines of QED in atomic units wqed = |(p- — 
, k +u2 ) ) the emission frequencies of the DFS in atomic units w DF = (137#) 3 + 137 3 # 2 (e 1 — e 2 ) 
and the values of the integer q of Eq. (40). 

In the DFS the interaction with a distant particle involves half the retarded Lienard- 
Wiechert potential plus half the advanced Lienard-Wiechert potential (henceforth called 
the semi-sum). This semi-sum yields a radiation magnetic field for the electron of (the 
far-magnetic field) 

(a_ x fj_) (a+ x n+) 

rad 2(l-n_-v_) 2 r 2(l + n+.v+)V 1 J 

where v and a are the electronic velocity and acceleration, n is a unit vector from the electron 
to the observation point, the subindex minus sign indicates evaluation on the retarded light- 
cone and the subindex plus sign indicates evaluation on the advanced light-cone. These 
two light-cones are defined by t± — t ± (r — n± • y), where y stands for the electron's 
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position. Along a precise circular orbit the first approximation to Eq. (47) has a zero 
spatial average. For the next term we avoid the Page expansion of Appendix A, because 
the deviating arguments are large; We approximate the size of Eq. (47) by expanding the 
denominators of Eq. (47), yielding the quadratic function 

B ,^ 2( fl .v)(ax fl ) (4g) 

We can estimate of Eq. (48) by noticing that along the n ± = x direction of the 
unperturbed plane this quadratic functional contains a product of the z perturbed coordinate 
times the x perturbed coordinate, i.e., the u and Z perturbations explained above Eq. (43). 
Translating the u mode to Cartesian coordinates with Eq. (36) we obtain 

B^oc^expW (49) 

According to Eq. (49), the frequency of the emission line is equal to Q plus the frequency 
of the uZ amplitude, 

w DF = Q + (e 1 -e 2 )Q/6, (50) 

with f2 given by Eq. (60). Notice that the emitted frequency of the DFS is naturally 
different from the orbital frequency. The fact that the emission frequency of hydrogen is 
different from the orbital frequency is a famous conundrum. The emission frequency of Eq. 
(50) contains differences of eigenvalues of the linear operator of Eqs. (39) and (34) and is 
strikingly similar to the Rydberg-Ritz combinatorial principle of quantum mechanics for the 
emission lines. 



VI. CONCLUSIONS AND DISCUSSION 



In the limit where the proton has an infinite mass, the concept of resonant dissipation 
loses meaning because the center-of-mass coordinate no longer plays a dynamical role. In this 
singular limit, there is a Lorentz frame where the proton rests at the origin at all times, and 
the field at the electron reduces to a simple Coulomb field in the DFS. The two-body dynam- 
ics in the DFS reduces then to the dynamical system of Eliezer's theorem; self-interaction 
plus a Coulomb field acting on the electron [3, 4]. We repeat this correct dynamics because 
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it is very unpopular [3-6] ; With inclusion of self-interaction, it is impossible for the electron 
to "spiral into the proton". Neither bound states nor dives are possible, only scattering 
states exist. This result is in surprising agreement with our formula (46) for the quantized 
angular momenta; If the mass of the proton is set infinite in Eq. (46), the quantized angu- 
lar momenta become infinite logarithmically, 9 goes to zero, and the particles are unbound 
at an infinite distance! One accomplishment of the present work is to recognize that only 
the two-body problem can produce a physically sensible electromagnetic-like model. Even 
though there is a dependence on the mass in Eq. (46), the logarithm of the mass ratio times 
6 A makes the theory very insensitive to this mass ratio, such that the deuterium and the 
muonium have essentially the same quantized angular momenta, in reasonable agreement 
with QED. Qualitative disagreement would need an exponentially massive charged particle. 
Fortunately to our present theory, such particle does not exist in nature. 

Another qualitative dynamical picture is suggested by Eliezer's result [3, 4]; The dynam- 
ical phenomenon that the electron always turns away from the proton along unidimensional 
orbits suggests that colinear orbits are the natural attractors of the dissipative dynam- 
ics (a ground state with zero angular momentum!). Along such orbits, the heavy particle 
(the proton) moves in a non-Coulombian way and the self-interaction provides the repul- 
sive mechanism that avoids the collision at the origin. This is again in agreement with 
the Schroedinger theory, where the ground state has a zero angular momentum. Again, 
the infinite-mass case produces unphysical dynamics; the electron turns away but then it 
runs away [4]. It remains to be researched if the two-body case has a physical orbit for 
zero-angular momentum orbits. 

The theory of normal forms for delay equations is studied in Ref. [27]. An analogous 
mathematical phenomenon is the finite-dimensional center-manifold for equations with ad- 
vance and delay studied in connection with discrete shocks in the conservation laws of Refs. 
[28, 29]. These conservation laws are similar to Dirac's relativistic Schroedinger's equation, 
and this would be a natural bridge to QED. Detailed construction of the resonant normal 
form is also needed to discuss the width of the emission lines. In the dynamical process of 
resonant exchange, the sharp line is emitted while the dynamics is locked to the neighbor- 
hood of the resonant orbit, which according to QED is a life-time of about 10 6 turns in the 
hydrogen atom (10~ 9 seconds). We conjecture that when the metastable orbit breaks down, 
the dynamics falls into the next metastable attracting orbit; another circular orbit, or into 
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the ground state [26]. 

The stiff modes of Eq. (40) introduce a fast (stiff) time-scale with a frequency of the 
order of a/6 ~ 1400 times the orbital frequency, such that the time for a stiff jump of the 
dynamics is times the orbital period, or 10 -18 seconds! After this fast timescale the 
resonance essentially locks the dynamics to the neighborhood of the metastable resonant 
orbit. The fact that the equations of electrodynamics describe stiff jumps in the phase 
space is largely unexplored in the light of modern applied mathematics; mainly due to the 
complexity involved. The dynamics starting from an asymptotic resonant orbit to another of 
a neighboring q is certainly described by a stiff jump, as expected generically from any stiff 
equation. In Ref [25], the much simpler Van der Pol oscillator is worked out in detail as an 
example of an equation of Lienard type that exhibits stiff jumps. In quantum mechanics one 
seems to need the problematic concept of an " instantaneous quantum jump" , to describe the 
stiff passage from one quantum state to another. It appears that classical electrodynamics 
prescribes exactly this qualitative phenomenon; a quasi- instantaneous fast dynamics. 

The dynamics in the DFS solves several conundrums of the classical hydrogen atom and 
is similar to QED in many ways; (i) the radiated frequency is not equal to the orbital 
frequency ( it is lesser than the orbital frequency, see Table 1). (ii) the resonant orbits are 
naturally quantized by integers and the radiated frequencies agree with the Bohr circular 
lines within one percent average deviation, (iii) The ratio of the emitted frequency to the 
orbital frequency is in reasonable agreement with QED. (iv) the angular momenta of the 
resonant orbits are naturally quantized with the correct Planck's constant, (v) the stability 
analysis uses a linear dynamical system with delay, a dynamical system that needs an initial 
function as the initial condition, just like Schroedinger's equation. The emitted frequencies 
are then given by a difference of two eigenvalues of this linear operator, like the Rydberg- 
Ritz combinatorial principle of quantum physics, (vi) The eigenvalues of our linear operator 
have a large magnitude that does not appear in the frequency. This large magnitude is given 
by a logarithm, just like in the divergent perturbation theory for the Lamb-shift of QED. 

Recognizing the correct qualitative dynamics with the concept of resonant dissipation has 
taken us very far; the stability analysis indicated the need for resonances, and these turned 
out to be satisfied only for the stiff modes and precisely in the atomic magnitude! The stiff 
modes also provide a natural integer to label the resonant orbits. We selected the values 
of q among the larger set predicted by the necessary condition (41), showing that Eq. (41) 
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is not in disagreement with QED. A sufficient condition should be part of the extra work 
to understand the unfolding of the bifurcation leading to the state of resonant dissipation. 
The large body of qualitative and quantitative agreement suggests that an extensive study 
of electromagnetic-like models [17], of which the DFS is only a generic example, could offer 
an explanation of QED in terms of a stiff dynamical system with third derivatives and delay. 
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VIII. APPENDIX A: PHYSICAL JUSTIFICATION OF THE DFS 

In Dirac's theory [1] the self-interaction is given by the sourceless combination of half of 
the retarded Lienard-Wiechert self-potential minus half of the advanced Lienard-Wiechert 
self-potential, i.e., the semi-difference [1]. This gives the following concise description of 
the DFS; Charges interact with themselves via the semi-difference of Lienard-Wiechert self- 
potentials and with other charges via the semi-sum of Lienard-Wiechert potentials. In the 
following we try to fit our ad-hoc DFS into Dirac's theory as an effect of the physical 
boundaries on the fields. Dirac's electrodynamics of point charges [1] uses the retarded 
potential F^ kret produced by each particle k and an incident free field F^ in . In Dirac's 
theory the electron and the proton of a hydrogen atom have the following equations of 
motion [1] 

2 2 

m 1 v lll --v llt --\\v 1 \\v llJi = -{F^ tin + F^ tret )v lv , (51) 
2 2 

m 2 v 2lM - -v^ - -\\v 2 \\ 2 v 2tl = (F£ in + F^ lret )v 2l/ , 

where double bars stand for the Minkowski scalar product, the electron and the proton have 
charges —1 and 1 respectively and the speed of light is c = 1. Since the DFS uses the 
semi-sum instead of the retarded-only potential, from the perspective of Dirac's theory this 
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demands the following constraints on the free field F^. 



F;, in (x 2 (t)) = \[p; 1M M{t)) - p; hret {x 2 {t))i 



(52) 



(53) 



where the field of each particle is to be evaluated along the trajectory of the other particle, as 
indicated by the parenthesis after each field. Since both the advanced and the retarded fields 
satisfy Maxwell's equations, the semi-difference is a free field, as assumed. The incident wave 
can be generated by the boundary conditions on the fields. For example, the reflections of the 
radiation by other atoms of a diluted gas could play the role of such a boundary condition. 

The semi-difference evaluated at the particle itself is the familiar self-interaction of the 
Dirac theory [1], and Eqs. (52) and (53) have instead the semi-difference evaluated at the 
position of the other particle. Using the Page expansion of the Lienard-Wiechert fields, we 
find that the electric field of this semi-difference is approximated by the third derivative of 
the other particle's coordinate, as discussed in Refs. [15, 16]. In this approximation with 
the Page series [15, 16], the incident electric field evaluated at the proton, Eq.(53), is 



Along the unperturbed orbit of Fig.l, Eq. (54) is an electromagnetic field rotating at the 
orbital frequency. For orbits in the atomic magnitude the electric field of Eq. (54) has an 
intensity that turns out to be of the order of the polarized vacuum of QED, as discussed in 
Ref. [31]. This shows that our needed homogeneous field has the correct physical magnitude 
of the QED vacuum polarized by the hydrogen atom. We see that the ad-hoc DFS demands 
a free field produced by the boundaries that is calculated to have a physically sensible order 
of magnitude. This approach to justify the DFS with a free field produced by the boundaries 
is similar to that of the stochastic electrodynamics of Refs. [32, 33]. 

Finally, we mention a more radical alternative to justify our ad- hoc DFS, by generalizing 
Dirac's theory such that the DFS would be derived from principle. This approach was taken 
by Eliezer and this generalization, henceforth called the Eliezer's setting (ES), is discussed 
in the excellent review of Ref. [17]. The ES involves the advanced interactions naturally, 
exactly in the same form of the DFS! Better still; the ES [17] contains an arbitrary parameter, 
and it would be highly desirable to experiment with stability analysis and the concept of 
resonant dissipation in the ES [17]. Even though the ES involves delay, advance and third 



E(x 2 (t)) 




2.. 



(54) 
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derivatives exactly like the DFS, the coefficients in the ES are never equal to those of the 
DFS. Our preliminary findings with the DFS suggest a future for this enterprise in the 
qualitative behavior of electromagnetic-like dynamics, one that could describe QED by a 
stiff dynamical system with delay. 



IX. APPENDIX B: DARWIN AND THE SOFT COULOMBIAN MODES 

In this appendix we calculate b\ and b 2 of Eq. (2) as a function of mi, m 2 and 9. The 
radial component of the electron's equation of motion along the circular orbit is [18] 

" >lWn,f) ~ = -j-{[C 2 + 9 2 S(S -!)](&! + & 2 cos(#)) + 9S(sm(9) - 9 cos(#))6 2 }, (55) 



where C and S are defined in Eqs. (13) and (14) respectively. Our Eq. (55) is Eq. (3.2) of 
Ref. [18] after use of Eq. (4) and the identity 

(1 - 9 2 b\)(l - 9 2 b\) =C 2 + 9 2 (S - 2)S. (56) 

The radial equation for the proton is obtained by exchanging the subindices 1 and 2 in Eq. 
(55). There are three equations involving bi, b 2 , 9 and r b ; (i) Eq. (55), (ii) the equation 
for the proton, obtained by exchanging indices 1 and 2 in Eq. (55), and (iii) the light-cone 
condition, Eq. (4) 

ii! j ° = tU[C 2 + 9 2 S(S - l)](&i + b 2 cos(0)) + 9S(sm(9) - 9 cos(9))b 2 }, 



: {[C* + z S(S- 1)](&2 + &i cos(9)) + 9S(sin(9)-9cos(9))b 1 }, 



m 2 b 2 r b 9 2 ■ , v ,., 2 



b\ + b\ + 2bib 2 cos(9) = 1. (57) 

For small values of 9 (atomic physics), we can solve Eqs. (57) in a power series of 9 with a 
symbolic manipulation software, yielding 

61 = w {1 + ^ )+d4D{mi ' m2) + - (58) 



where 



n/ \ (ft xr 12m i _ 13m i - 5mim| + \\m 2 m\ 

D( mi ,m 2 ) ee (_)[ — ]. (59) 
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It is easy to continue this power series, but for the stiff limit in the atomic magnitude, even 
the 8 2 correction already gives a very small correction. The orbital frequency is determined 
by 

the first term is Kepler's third law if we use 6 = Qr&, and the next term is the Darwin 
correction. More information about the isolated two-body problem can be found in Refs. 
[21, 23]. 

As an application of the above expansion, we calculate the soft Coulombian modes of Eq. 
(31) at a finite p by expanding up to 0(9 3 ) : 

M 9xr 6 2 , 12/i sl 2M o, , a . . . 

det Z = -p\l + p 2 )[l - _(l - J)] - g— pV + 2^). (61) 

The soft roots of det Z = for Eq. (61) are the Galilean translation mode p = (a 
double root) and the oscillatory solutions p ~ ±i that have a real part describing the 
radiative damping of the DFS, a familiar feature. We had partial success describing the 
atomic dynamics of helium with the Darwin approximation [30], and the tools of stability 
analysis used here were already used in Refs. [15, 16]. The concept of resonant dissipation 
is new, and it is a generalization of the concept of a non-ionizing dynamics of Ref.[30]. 
Unfortunately, the theory of Refs. [15, 16, 30] fails to describe discrete states for hydrogen 
because the soft Darwin modes are neutrally stable. As we have seen here, it is the stiff 
modes that equilibrate the dynamics, and those are beyond the Darwin approximation. 

X. CAPTIONS 

Fig. 1: The unperturbed circular orbit with the particles in diametral opposition at the 
same time in the inertial frame. Indicated is also the advanced position of particle 2 and the 
angle travelled during the light-cone time. The drawing is not on scale; The circular orbit 
of the proton has an exaggerated radius for illustrative purposes. Arbitrary units. 
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